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Abstract 

A simple and intuitive calculation, based on the semiclassical approximation, 
demonstrates how the large size of a hadronic target introduces a new pertur- 
bative scale into the process of small- x deep inelastic scattering. The above calcu- 
lation, which is performed in the target rest frame, is compared to the McLerran- 
Venugopalan model for scattering off large nuclei, which has first highlighted this 
effect in the infinite momentum frame. It is shown that the two approaches, i.e., 
the rest frame based semiclassical calculation and the infinite momentum frame 
based McLerran-Venugopalan approach are quantitatively consistent. 



1 Introduction 



It has been suggested by McLerran and Venugopalan |T[ that for very large nuclei the 
parton distribution functions at small x are perturbatively computable. Their argument 
relies on the high density of partons per unit area per unit rapidity which serves as a 
new hard scale in the problem. Extending these ideas, an evolution equation has been 
formulated || to follow the growth of this potentially large density as x decreases. These 
calculations are formulated in the infinite momentum frame of the hadron. 

In an apparently quite different approach, small-x deep inelastic scattering is de- 
scribed in the target rest frame, where the long-lived partonic fluctuations of the very 
energetic photon scatter off the target hadron ||. Explicit formulae focusing on the qq 
component of the photon can be found in Q]. In the semiclassical framework, the scatter- 
ing of the partonic fluctuations off the target is modeled by an eikonal interaction with 
the target colour field |5|-|7] . 

In the present paper, the semiclassical approach in the target rest frame is applied to 
the case of very large hadronic targets, such as heavy nuclei. The perturbative calculabil- 
ity of parton densities advertised previously emerges in a very simple and intuitive way. 
It is the result of the limited transverse size of the relevant qq fluctuations of the virtual 
photon. For an ordinary hadronic target, leading twist contributions to deep inelastic 
scattering come from both small and large size qq configurations. For larger targets, the 
qq pair has to penetrate an extended colour field and the geometric limit of the cross 
section is already reached for small size pairs. Accordingly, non-perturbative large size qq 
configurations do not contribute significantly to the total photon hadron cross section. 

We then go on to compare the semiclassical approach with the McLerran- Venugopa- 
lan approach in greater detail. It is shown that they can be mapped onto each other once 
the non-Abelian eikonal factors characteristic of the target colour field are recognised 
as the basic building blocks of both methods. We confront the different perspectives 
inherent to the rest frame and infinite momentum frame calculation. 

Performing an average over different colour field configurations, explicit formulae are 
obtained for the dipole cross section a(g), which describes the interaction of a qq pair of 
transverse size g with the target. 

Although the above results are only valid in the limit of very large targets, we believe 
they are useful in a much wider sense. They provide a theoretical laboratory in which, 
starting from a non-perturbative yet technically manageable description of the target, 
different small-x processes can be investigated. In particular, they establish a non-trivial 
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limit in which the semiclassical approach can be justified in the framework of QCD. 

The paper is organized as follows. After summarising the target rest frame view on 
deep inelastic scattering, a simple derivation of the perturbative scale associated with 
large targets is given in Sect. Q The equivalence to the McLerran-Venugopalan approach 
is demonstrated in Sect. |^ and explicit models for the averaging over the gluon field are 
described in Sect. [|. The conclusions are given in Sect. ^. 

2 Calculation in the target rest frame 

The total cross section for the scattering of a virtual photon off a hadronic target can be 
obtained from the Compton scattering amplitude via the optical theorem. In the target 
rest frame, the leading order partonic process is the fluctuation of the photon into a qq 
pair (see Fig. [[]). To keep calculations as simple as possible, we consider a longitudinally 
polarized photon coupled to scalar quarks with one unit of electric charge. As far as the 
Q 2 -behaviour is concerned (Q 2 = —q 2 ), this is analogous to the standard partonic process 
where a transverse photon couples to spinor quarks f|. 




target 



Figure 1: The Compton scattering amplitude within the semiclassical approach. 

In analogy to @], the longitudinal cross section can be written as 

a L = J d 2 g L a(g)W L (g) , (1) 

with the square of the wave function of the virtual photon given by 

W L (g) = ^JdaN 2 K 2 (Ng). (2) 

Here g = \g±\is the transverse size of the qq pair, a is the longitudinal momentum fraction 
of the photon carried by the antiquark, N 2 = a(l — <y)Q 2 , and Kq is a modified Bessel 
function. Note that, in contrast to J3|, we have defined Wi to include the integration 
over a. 

The cross section for the realistic case of a transverse photon and spinor quarks is 
obtained by substituting K 2 (gN) with 2[a 2 + (1 - a) 2 ]K 2 (gN) in Eq. (|). None of the 
qualitative results derived below is affected by this substitution. 
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Within the semiclassical approach, the dipole cross section a(g) is given by 



a(g) = — j d 2 x±tr 



Q± 



(3) 



where the SU(3) matrices U and W represent the non-Abelian eikonal factors associated 
with the quark and antiquark propagating through the target colour field and averaging 
over all field configurations is implicit. 

The functional form of cr(g) is shown qualitatively in Fig. |2]. For conventional hadrons 
of size ~ 1/A (where A ~ Aqcd) its typical features are the quadratic rise at small g 
(cr(g) ~ g 2 with a proportionality constant 0(1)) and the saturation at a(g) ~ 1/A 2 , 
which occurs at g ~ 1/A. Consider now the idealized case of a very large target of size 
77/ A with i] ^> 1 (77 ~ A 1 / 3 for a nucleus). It is easy to see that at small g the functional 
behaviour is given by a(g) ~ rfg 2 while saturation has to occur at a(g) ~ i] 2 /A 2 for 
geometrical reasons. It follows that the change from quadratic rise to constant behaviour 
takes place at g ~ l/y^A, i.e., at smaller g than for conventional targets. 




large hadron 
(radius ~r|/A) 



conventional 
hadron 
(radius -1/A) 



> e 



-l/frjA 



Figure 2: Qualitative behaviour of the function a(g). 

From the above behaviour of a(g) we will now derive the dominance of small trans- 
verse distances in the convolution integral of Eq. (H). To do this, a better understanding 
of the function Wi{g) is necessary. Recalling that Kq(x) ~ ln(l/x) for i«l while being 
exponentially suppressed for x 3> 1, it is easy to see that Wi{g) ~ Q 2 ln 2 (l/^ 2 Q 2 ) for 
g <^\/Q and Wi(g) ~ l/f? 4 Q 2 for g ^> 1/Q. Here numerical constants and non-leading 
terms have been suppressed. 

Under the assumption A 2 77 A 2 <ti Q 2 , the integral in Eq. (JJ) can now be esti- 
mated by decomposing it into three regions with qualitatively different behaviour of the 
functions Wi{g) and cr(g), 

VL = al + 4 + v?=( +/ /2 +/ Udg 2 a(g)W L (g). (4) 
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Of the three contributions 



/•VQ 2 n 3 

I / J„2„3„2 s^2i„2/i l n 2r\2\ 'J_ 

Q 2 



f d Q 2 rf Q 2 Q 2 \n 2 {l/ Q 2 Q' 
Jo 



ii/ V A' a 2 ^ 4 g 2 g 2 

the second one dominates, giving the total cross section 

^~|>(Q 2 M 2 )- (6) 

It is crucial that the third integral is dominated by contributions from its lower limit. 
Therefore, the overall result is not sensitive to values of g that are larger than 1/^/r/A. 
Phrased differently, for sufficiently large targets the transverse size of the qq component 
of the photon wave function stays perturbative. 

3 Calculation in the infinite momentum frame 

The above considerations can be related to the infinite momentum frame calculation 
of parton distributions within the McLerran-Venugopalan model ||. Here we use the 
term infinite momentum frame for a reference frame where the '+'-component of the 
hadron momentum and the '-'-component of the photon momentum are large (k± = 
(k ± fe 3 )/v2). We start with the operator definition of scalar parton distributions 

9(0 = — [dy.dx.d 2 x ± e-^ p+y - <P|T0 t (O, x_ + y_, x ± )U x _+y_, x _<j)(0, x_, x ± )\P> . 

7T J 

(7) 

Note, that we use a time ordered product in this definition (see e.g. flOfl ) and choose to 
normalize the hadron state by demanding < P\P >= 1. The matrix U is the standard 
link operator required for gauge invariance of the above definition. The corresponding 
space time picture for a hadron localized at z_ = is illustrated in Fig. |3|. 

It is understood that only the connected diagrams contribute to the definition of g(£) 
in Eq. (0), which is equivalent to subtracting the vacuum expectation value of It 
follows that we just need to calculate the propagator of the scalar field in the presence 
of the hadron and to subtract the propagator in the vacuum, 

q(£) = ^—t[ dy.dx^d 2 x ± e- iiP+y -ti [U(x ± )G(x, x + y) - G (x, x + y)} . (8) 

71 J 
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Figure 3: Space-time picture in the infinite momentum frame. 

Here x = (0, x_, x±), y = (0, y_, 0) and Go is the conventional Feynman propagator. 

At very small x, the contribution of the constituent quarks to the quark distribution 
is negligible. It is therefore sufficient to calculate the propagator in the background of a 
certain colour field and to average over all field configurations of the hadronic state. In a 
gauge in which the potential vanishes outside the hadron (say, outside the shaded region 
of Fig. |H) the propagating particle is simply rotated in colour space at the moment when 
it penetrates the field. The high energy limit makes the region with non-zero potential 
arbitrarily narrow in z_-direction so that, for ?/_ > and X- < 0, the propagator can be 
written as 



G(x,y) = Jd 4 zGo(x,z) 



25(z„)U\z L 



d 

idz + 



G (z,y). (9) 



Note that this corresponds to the propagation of a positive energy antiquark from x to 
y. The non-Abelian eikonal factor at transverse position z± is explicitly given by 

/+oo 
dz-A+(z-, z±) , (10) 
-oo 

where is the gauge potential of the hadron. An expression analogous to Eq. (Bl) holds 
in the region y_ < and X- > 0, while for the remaining regions where x_ and y_ 
are either both positive or both negative the propagator is simply the free one. This 
intuitive physical picture is slightly less apparent if one uses an A + = gauge as in ||. 
Then the propagators contain an additional eikonal factor corresponding to the relative 
gauge transformation while the link operator in Eq. (|7|) is absent. 

The expression for the propagator provided by Eq. (pj) is exact for a 5-function-like 
source. Within the present context, the highly Lorentz contracted target has a 5-function- 
like appearance on scales larger than 1/P+. The use of Eq. ([|) for the quark density 
calculation according to Eq. (|J) is justified since typical values of ?/_ ~ l/£-f + are much 
larger than 1/P+. 
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Note furthermore, that the gauge potential in Eq. (|T0D is extremely large in the 
infinite momentum frame. However, this is compensated by the smallness of the relevant 
integration region in z_. Clearly, the explicit formula for the eikonal factor Eq. ([H]) is 
invariant with respect to arbitrary boosts along the Z3-axis. 

Inserting our expression for the propagator into Eq. and performing the necessary 
integrations || the following result is obtained [] 

= 7A4 NW J d 2 ^K 2 (Ng)a(g) . (11) 

The dipole cross section a(g) (see Eq. @) emerges naturally due to the presence of two 
eikonal factors from the combination of Eqs. (|7]) and (§). The ultraviolet cutoff Auv < ^ 
Q is introduced to suppress the contribution from partons with very high transverse 
momenta (compare the discussion in [llj). 

Next, we want to explain in which sense Eq. (O) is equivalent to the cross section 
formula Eq. ([I]). The cross section <jl receives leading twist contributions from two re- 
gions, the aligned jet region, where a <C 1 or (1 — a) <C 1 and g is large, and the high-p^ 
region, where a ~ 1/2 and g is small. The small-a contribution of Eq. ([I]), 

* L * = ^ ^ I d2 ^ N2K o( N ^) • ( 12 ) 

corresponds, in parton model language, to the contribution from the quark distribution 
of the target (hence the index q). The latter can be identified by comparing with the 
parton model formula (for scalar partons) 

o L , q = g 2 xq{x) . (13) 

As expected, the obtained quark distribution is in agreement with Eq. (pTTj) . 

The perturbative nature of the scattering process off large nuclei, as derived in Sect. 
from generic features of cr(g), is implemented within the McLerran-Venugopalan model 
via the colour field average implicit in Eq. (|IID . This is subject of the next section. 

4 Explicit models for a(g) 

Consider first a particularly simple model that can be analysed with minimal calcula- 
tional effort. 

1 Note that a sign error in Eq. (47) of M leads to a spurious cancellation which affects their result. 
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According to Eq. ([J), the expectation value of UW in the background of a hadron is 
required. The averaging over the field configurations of an ordinary hadron with radius 
~ 1/A is symbolized by (• - -}i/a- For values of x± that correspond to a central collision 
and in the limit of small g, we have 

( UaixMixj. + g±) ) 1/A ^ M« " Bq 2 K 2 T^ . (14) 

Here k, I are colour indices, the matrices T a form the conventional basis of the SU(3) 
Lie algebra and B is a constant. The simple colour structure follows from the fact that, 
at small g, the path ordered exponentials implicit in UU> can be expanded in powers of 
the gauge field A a . The leading contribution comes from (A a A b ) ~ 5 ab . 

To extend Eq. flU) to the case of a large hadron, an additional assumption has 
to be made. It is assumed that the volume of the hadron can be split into mutually 
uncorrelated regions of size ~ 1/A, as the naive geometrical picture of a large nucleus 
built from nucleons suggests. If the linear extension of the hadron is ~ 77/A and n ~ rj 
is an integer, the qq pair passes approximately n such regions on its way through the 
hadron (see Fig. [I]). The expectation value of UW is given by 

( U{x ± )rf(x ± + g ± ) ) v/A * ( ( U(x ± )rf(x x + g ± ) ) 1/A J • (15) 

Here the appropriate contraction of the colour indices on the r.h. side is implied. Using 
Eq. flTJD the trace of Eq. ([14]) can be calculated for the large hadron case. For sufficiently 
large 77 and small g, the product of C/'s exponentiates giving 



tr 



U(x ± )U j (x ± + g ± )] ) A ~N c (l-C F Bg 2 A 2 y~N c exp[-r]C F Bg 2 A 2 } (16) 



n/A 

with Cf = (N% — 1)/2N C . Again, x± is chosen such as to describe a central collision. 
Therefore, the above is valid for an x± area of size ~ r] 2 /A 2 . Inserting this into Eq. (|3]) 
the following result is obtained 

cr(g) ~ j^N c (l - exp[- V C F Bg 2 A 2 }) . (17) 

Clearly, this simple model is consistent with the generic features of a(g) discussed in 
Sect. H The above expansion in powers of the gauge potential is only valid if g <C 1/A. 
The exponentiation in Eq. (|16|) is valid in the limit of large rj for g ~ 1/y/rjA. Both 
conditions are fulfilled in deep inelastic scattering in the limit of very large hadronic 
targets. The simple formula of Eq. (|i~7"D, which includes the numerical constant B and 
an unknown overall normalization, provides a model for a(g) in this limit. 

Let us now consider the same quantity within the framework of [|l],|]]. Here, one of the 
central issues was to provide an explicit prescription for the averaging over the target's 
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colour fields. The initial suggestion, the McLerran-Venugopalan model, used a Gaussian 
statistical weight for the colour charge in the target that generates these fields (compare 
also the discussion in |I2[). Later, an evolution equation was derived to determine such 
a statistical weight and its x dependence 0. These results suggest that, although a 
Gaussian weight is not a solution of the evolution equation, it may be used at not too 
small x, as long as the parton densities are not too high. In the present calculation we 
satisfy ourselves with the simple Gaussian weight introduced in 

As in [0,0, we calculate the classical field for a static, i.e. x + independent, colour 
source of the form 

J«(x)=p a (x_,x ± )g,_, (18) 

where p parametrises the adjoint charge density of the hadron (not to be confused with 
the size g of the qq pair used above). In the A_ = gauge the corresponding gauge field 
is given explicitly by 

V1A + (x_,x_l) = gp(x_,xj_) ; A ± (x-, x±_) = , (19) 

where g is the gauge coupling. 

Note, that here we are concerned with the space time region where the potential in 
Eq. fllCf ) is nonzero. From the point of view of Eq. @, where much larger distances in 
x„ dominate, this corresponds to resolving the structure of the ^-function in x_. 

According to |IJ , the Gaussian average over p is implemented independently at each 
position in x_. This means that the average of an arbitrary functional 0[p] is given by 




where the parameter p? ~ A 3 determines the typical value of the colour charge density. 
The above average implies a correlation function 

(A°.(x-,x ± )A%(y-,y ± )) = g 2 p 2 6 ab 6(x- - y_) 7 (*L - y±) (21) 

within the target. Here 7 is the inverse of the differential operator in Eq. (p0| ) and is 
given by 

1 . . r d 2 k ± e ~ ik ^ . . 

7( ^ ) = vi (x±) = A(^F75F' (22) 
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where an infrared cutoff A has been introduced to define the k±_ integral. It follows that, 
at sufficiently small x±, 

7(^)-7(0)^hi(xiA 2 ). (23) 

Non-logarithmic terms depend on the specific way of introducing the cutoff and have 
been neglected in the above equation. 

Using Eqs. ( p0| ) and (|2T|), the required average of a product of eikonal factors can be 
evaluated along the same lines as the gluon distribution in JIB] , 

(tr [U(x ± )tf(x ± + Q±)]) v/A = N c exp (g 4 C F fl 2 [ 7 (f?x) - 7 (0)]) (24) 

where p 2 = J dx^p 2 ~ rjA 2 reflects the large linear extension 77 of the target. According 

target 

to Eq. (|3]), an x± integration has to be performed over a transverse area of a size ~ 77 2 /A 2 , 
giving rise to the formula 



,2 



77 

o(q) ~ ^2 Nc ( 1 ~ exp 



-^^A 2 , 2 ln(l/, 2 A 2 ) 



(25) 



Due to the simple model for 7 used above, this formula is only valid for p <C 1/A. 
However, as has been discussed in Sect. |2|, this is sufficient to determine the virtual 
photon cross section for large rj. The result shows the same Glauber type (1 — exp) 
structure as Eq. ( |TT| ) and the gluon distributions calculated in [ |T3"| , [T4|j . 

The present analysis is concerned with an energetic colour dipole created far way 
from the target and scattering off the target colour fields. It is the combination of the 
propagation of this dipole, as encoded in the Bessel functions in Eq.(g) or the Feynman 
propagators of Eq. (^), and its interaction with the target that determines the quark 
distribution. A similar analysis could be performed replacing the pair of quarks with 
a pair of gluons. This would result in a contribution to the gluon density beyond the 
instantaneous interactions considered in [O . 



5 Conclusions 



The scattering of a virtual photon off a large hadronic target, such as a very large nucleus, 
has been considered in the high energy limit. The process has been discussed using both 
the semiclassical target rest frame calculation and the McLerran-Venugopalan approach 
in the infinite momentum frame. Both approaches show that, for sufficiently large targets, 
the cross section is dominated by perturbative contributions. 

More specifically, the cross section can be described in terms of the qq component 
of the photon scattering off the target. For large targets, the process receives no lead- 
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ing order contributions from quark pairs with large transverse size. This means that 
higher Fock states of the photon are expected to be suppressed by powers of a s and the 
semiclassical approach can be justified in QCD. 

One of the main problems of the presented approach is certainly the procedure of 
determining the target colour fields which form the non-perturbative input of the calcu- 
lation. Two simple models for the averaging over these fields have been outlined but no 
satisfactory QCD based procedure is presently known. 

We would like to thank W. Buchmiiller, A. Kovner, L. McLerran and A.H. Mueller for 
helpful discussions and comments. 
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